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In one-dimensional case the search for presence of the anomalous phenomena in 
multiplicity distributions is usually performed in frame of the horizontal, vertical 
and mixed types of the analysis. We show that if the data involve a d-dimensional 
phase space, there exists a convenient procedure, generalizing the one-dimensional 
approach, which allows to introduce more non-trivial types of the analysis; we 
formulate them in some detail in one-, two-, and partially also in three-dimensions. 



1 Introduction 

Search for anomalous phenomena at high energies is usually performed in terms 
of a convenient kind of statistical moments, the information about the presence 
of those phenomena being encoded in appropriate sort of the scaling indices. 
Quite often, the one-dimensional cases are considered 1 and usually three types 
of the analysis are applied there 2 , namely, the horizontal (H), vertical (V), and 
the mixed (HV) ones. However, there is still missed a systematic approach to 
higher dimensions as well as a deeper rooted physical theory or at least a more 
or less reliable dynamical model for the scaling indices characterizing those 
phenomena. 

Let us add that in the cases considered in the present contribution, the 
appearance of the non-statistical, anomalous fluctuations, induces the need to 
take into account the non-continuous statistical distributions. 

In principle, the notion "anomalous" means (special, see below) "non- 
statistical" while "statistical" denotes "pure Poissonian" . 

When investigating the presence and properties of the anomalous phenom- 
ena in high energy physics, the factorial moments of rank q are usually applied, 
their normalization being such that for the pure Poissonian distribution they 
are equal to unity (for q= 1, 2, 3, ...). This means that their deviation of 
unity (at least for some values of the rank q) signifies the presence of the 
non-statistical behavior in the corresponding data. 

In the present contribution a sample of E events is investigated, the count- 
ing index of the events being denoted by e, i.e. e=l, 2, 3, E. The number 
of (charged) particles observed in a given region of the phase-space under con- 
sideration is denoted by n. In what follows also the factorial multinomial, T, 
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introduced by 



T = T{n-q)= n(n-\)...(n-q+l) , (1) 

is applied. 

2 One-dimensional case 

We assume that the dependence of the multiplicity distribution of (charged, 
produced) secondaries on one variable (say, the pseudorapidity 77 of every par- 
ticle under consideration) is known from the observed collisions. 

Let the multiplicity data be known in a pseudorapidity window A77, and, 
eventually, this window be partitioned into M equal size bins, 5t], Si] = (Arj) /M; 
the index counting the bins is denoted by m, i.e. m = 1,2,3,..., M. It is un- 
derstood that the number of particles in every elementary cell specified by two 
numbers, namely, (me), i.e. n mei is deduced from the experimental output. 

Usually it is said that the anomalous phenomenon (intermittency) is present 
if the experimental data exhibit (with increasing number of bins, M) a linear 
dependence between the logarithm of the factorial moments and logarithm of 
the number of bins, M, i.e., if the factorial moments reveal a singularity of 
the form M a " , a q being non- vanishing quantities; they are called slopes or 
scaling indices 3 . It is worthwile to introduce also the following notations: to 
emphasize the Individual, Independent averaging of the expression entering 
the numerator or denominator of the corresponding moment, 

1 M 1 E 

Im e mE ' ll = e E ' 

m— 1 e— 1 

and, to emphasize the averaging over Both, the numerator as well as the 
denominator, 

1 M 1 E 

m— 1 e— 1 

Now, application of the horizontal type of the analysis means that the fac- 
torial moments in the (me)— th cell [F(q; M)]me , m the e— th event [F(q; M)]i H ^ 
and in the whole sample [F(q;M)]^ H \ are introduced, respectively, by 



(2) 



and 

[F(«;M)]W = I m [F(q;M)}W 
_ l m F{n = n me ;q) 



(3) 



[I m ' Tlm'e]^ 

[F(q;M)]W = B e [F(g;M)]W (4) 

where J 7 is given by Eq. (1). In this case the scaling properties are investigated 
by means of the relation, 

[F{q;M)]W cx f^M< H) . (5) 

The vertical type of the analysis is specified by the following form of the 
factorial moments, 

|± e ' ll m e'\ 

and 

[F{q-M)]<y) = I e [F(q;M)]<y) (7) 

[F(q;M)]W = B m [F(q;M)]W 

= BMF(q;M)]W (8) 

where T is still given by (1). Now, the scaling property is formulated in the 
following way, 

[F{q;M)]W cx /f>M<> . (9) 
The mixed type of the analysis is characterized by the following expressions, 

[±m' *-e' li "m'e' \ 

and 

[F(q;M)r^ = I m I e [F(q; M)}^ oc f^M^ . (11) 
We conclude this section by the observation that the expressions 

B e I m , B 7n I e , I m I e 

represent, respectively, a shorthand notation of the three types of the analysis 
mentioned above. It is perhaps worth to mention that in case of the silicon (at 
14.6 A GeV/c) colliding with the emulsion nuclei (the BNL data), those three 
types of the analysis lead 4 to similar values of the slopes (scaling indices) a q , 
with q = 2,3,4,5. 



3 Two-dimensional case 



Let us investigate the presence of non-statistical fluctuations in the (charged) 
multiplicity distributions depending on two variables (like, e.g., the pseudo- 
rapidity and azimuthal angle). We present shortly the way which allows to 
formulate the relations appropriate to this case. 

First of all, we introduce the partitioning of the windows where both quan- 
tities mentioned above were measured. The corresponding bins are counted by 
the indices m\ and m 2 where 

1 < mi < Mi and 1 < m 2 < M 2 . 

The number of charged particles, n, observed in the (mi, m 2 , e)-th elementary 
cell is denoted by 

^ — ^mi,m2,e • (1^) 

In this case the following notation is useful, 

Imi = II j Im 2 = I2 ! B mi = Bi , B m2 = B 2 

and I e , B e represent the same averaging as in the preceding Section. 

There are seven non-trivial types of the analysis; they can be represented 
in the following form, 



[2dim; 1] 
[2dim; 2] 
[2dim; 3] 



B e B 2 Ii, [2dim;4] 

B e B!l 2 , [2dim;5] 

B 2 BiI e , [2dim;6] 

[2dim; 7] : I e I 2 Ii . 



B e I 2 Ii 
B 2 I e Ii 

B!l e I 2 



For instance, in the first case, [2dim; 1], 



T{n\ q) 



[Im]i 



Iif(n;q) 

[lm]i 



F(q) = B e B 2 ^^ cx .fq(MiM 2 ) aq 

where the multiplicity n of the elementary cells is given by Eq. (12). 
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(13) 



(14) 



4 Three-dimensional case 



If the question is to be answered whether the non-statistical fluctuations are 
present in multiplicity data depending simultaneously on three variables (like 
e.g. the pseudorapidity, azimuthal angle and transverse momentum), first of all 
the number of (charged) particles in every elementary cell should be specified. 
To this end, the windows in all three variables are partitioned thereby involving 
the relations, 1 < nij < Mj 7 with j = 1, 2, 3. In this case, the multiplicity n of 
the elementary cell is specified by 



n = 



i, mi,m2,m3,e 



(15) 



Introducing the notation analogous to that one applied in the preceding cases, 
namely 

Imj = Ij an d B mj = Bj j = 1,2,3, 

we observe that in the present case the following fifteen possible types of the 
analysis can be introduced, 



[3dim;l] : B e B 3 B 2 Ii 
[3dim;2] : B e B 3 BiI 2 

[3dim;5] : B e B 3 I 2 Ii 
[3dim;6] : B e B 2 I 3 I 1 
[3dim; 7] : B 3 B 2 IeIi 

[3dim; 11] : B e I 3 I 2 Ii 
[3dim; 12] : Bglefcli 



[3dim;3] : B e B 2 B 1 I 3 
[3dim;4] : B 3 B 2 B!l e 

[3dim;8] : B e BiI 3 I 2 
[3dim;9] : B 3 BiI e I 2 
[3dim; 10] : B 2 BiI e I 3 

[3dim; 13] : B 2 I e I 3 Ii 
[3dim; 14] : B^Igl, 



(16) 



[3dim; 15] : I e I 3 I 2 Ii . 

In Eq. (16), the B e and I e have the same meaning as in the preceding Sections. 
For instance, the first type of the analysis mentioned in Eq. (16) involves the 
factorial moments in the (mi, to 2 , m 3 , e)-th cell, 

_ T{n-q) 

as well as the full form of that moment together with the corresponding scaling 
condition, 

F(g)=BeB 3 B2liF roi , mai m3,e(«) « f q (M 1 M 2 M 3 )^ (18) 

5 



where the multiplicity n is given by Eq. (15). 

Let us add that every type of the analysis is characterized by its own set of 
the scaling indices a q and the corresponding form of all kinds of the associated 
moments (like the frequency, correlation, dispersion, etc. moments); some 
more details in 5 . 



5 Dispersion moments 

It is well known that several multiplicity distributions observed in the past, 
were satisfactorily described by continuous statistical distributions inherent 
in different modelling ideas. And quite often the associated moments, espe- 
cially the corresponding dispersion moments, represented a very apt tool which 
helped to discriminate between more and less appropriate models. Now, it can 
be expected that especially at higher energies (available today as well as in 
the future) the possible presence of the intermittency should be taken into ac- 
count. In this case, application of the continuous distributions is already not 
sufficiently adequate, compare in particular e.g. Fig. 1 in 6 ; that distribution, 
in 6 , representing one event, is induced by non-self-similar processes and clearly 
reveals scaling property 7 (also other authors deal with the individual events, 
for instance 8 ). Application of the corresponding associated statistical moments 
(like e.g. the appropriately introduced dispersion moments), might play again 
a decisive role in the procedure leading to the discern between the more or less 
reliable fractal model approaches. Moreover, the dispersion moments involve 
also a piece of information about the presence of the empty bins. 

In particular, the first type of the analysis mentioned in the present paper, 
[ldim; 1] , involves (let us call them "modified" ) dispersion moments D defined 
by the following way 

[*e=(j^-l) q (19) 

and D q = B e I m [D|] me . On the other hand, the "standard" dispersion mo- 
ments D are introduced by 

[Dq}me ij^rne Im'^m'e)^ i (20) 

i.e. their relation with the modified D-dispcrsion moments acquires the form 

= , 1 ]q [JgU (21) 

and 

D q q = B e Im[ ^ ]m ; = BeM^U • ( 22 ) 

|± m ' ii>rn'e\ 
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Moreover, in the case [2dim; 1]: 



Pj]mi,n»2,e — ( T n I ' ^9 B e B m2 I mi [D^] mi m2 e ; (23) 

and [3dim; 1]: 



l' 'H"m' 1 ,m2.e 



rrwi I '"mi ,m,2 ,m,3 ,e _ , 



Jg=B e B m3 B m2 I mi | _ -mi,m 2 ,m 3 ,e _^ (24) 

We note that the multifractal analysis (of the type [ldim; 1]) of data 
coming from the collisions of gold (at 10.6 A GeV/c) with emulsion nuclei 9 
leads to the conclusion that the lower rank dispersion moments (essentially of 
the form of Eq. (19)) reveal scaling property with sufficient accuracy 10 . 

Analogous procedure can be applied also in higher dimensional cases. 



6 Conclusions 

In the present contribution, investigating presence of the anomalous fluctua- 
tions, a natural and systematic extension of the one-dimensional approach to 
arbitrary, say, d-dimensional case is proposed; it involves 




= 2 d+1 - 1 



non-trivial types of the analysis. The aim is to apply such a systematic ap- 
proach to a same set of the data and to look for the regularities embedded in 
the scaling indices which belong to the same branch of the analyses. Hope- 
fully, the regularities which would be observed there might lead at least to the 
construction of the corresponding model equations. 

Especially, it is expected that a systematic application of the non-trivial 
types of the analysis to several samples of events will allow (i) to recognize 
which type of the analysis is the most suitable one for detecting the presence 
of the concrete type of the anomalous fluctuations, (ii) to introduce an appro- 
priate cathegorization of the anomalous phenomena, (iii) to ascribe a concrete 
numerical value to the parameters which enter theoretical expressions applied 
to the description of the experimental data, and, eventually, (iv) to formu- 
late a sufficiently reliable (and still missing) theory (or at least model) of the 
anomalous phenomena. 
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